Constraining quark angular momentum through semi-inclusive measurements 
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The determination of quark angular momentum requires the knowledge of the generalized par- 
ton distribution E in the forward limit. We assume a connection between this function and the 
Sivers transverse-momentum distribution, based on model calculations and theoretical considera- 
tions. Using this assumption, we show that it is possible to fit at the same time nucleon magnetic 
moments and semi-inclusive single-spin asymmetries. This imposes additional constraints on the 
Sivers function and opens a plausible way to quantifying quark angular momentum. 
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Nucleons are spin-1/2 composite particles made by 
partons (i.e., quarks and gluons). Determining how much 
of the nucleons' spin is carried by each parton is a crit- 
ical endeavour towards an understanding of the micro- 
scopic structure of matter. In this work, we propose a 
way to constrain the longitudinal angular momentum J° 
of a (anti)quark with flavor a. To do this, we adopt an 
assumption, motivated by model calculations and the- 
oretical considerations, that connects J" to the Sivers 
transverse-momentum distribution (TMD) measured in 
semi-inclusive deep-inelastic scattering (SIDIS) ^1 . The 
Sivers function /-j^ [5] is related to the distortion of the 
momentum distribution of an unpolarized parton a when 
the parent nucleon is transversely polarized. We show 
that this assumption of relating J° to f^ is compatible 
with existing data, and we derive estimates of J"" . 

The total longitudinal angular momentum of a parton 
a (with a — q,q) at some scale Q^ can be computed 
as a specific moment of generalized parton distribution 
functions (GPD) [31 

J'^(q2)^1 f dxx(H''{x,0,0;Q'^)+E''{x,0,0;Q'^) 
2 Jo ^ 

(1) 
The GPD H''{x,0,0;Q^) corresponds to the familiar 
coUinear parton distribution function (PDF) fi{x;Q^), 
which gives the probability of finding at the scale Q^ 
a parton with flavor a and fraction x of the (longitu- 
dinal) momentum of the parent nucleon. The forward 
limit of the GPD i?° does not correspond to any coUinear 
PDF [4 . It is possible to probe the function E'^ in 
experiments, but never in the forward limit (see, e.g., 
[5]). Assumptions are eventually necessary to constrain 
£'"(x, 0,0; Q2). This makes the estimate of J° partic- 
ularly challenging. The only model-independent con- 
straint is the scale-independent sum rule 
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dxE^-^x, 0,0) = K, 



(2) 



where E"^^ = E'^ — E^ and k denotes the anomalous mag- 
netic moment of the parent nucleon. 

Inspired by results of spectator models [BHTU] and theo- 
retical considerations [T] , we propose the following simple 
relation at a specific scale Ql, 

ft'^^^^x; Ql) = ~L{x) E-{x, 0, 0; Ql), (3) 

where we define the n-th moment of a TMD with respect 
to its transverse momentum k± as 



ft^^^^^{x;Q')^ Id'k^l^^ 



ftT^{x,ki;Q'), (4) 
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and M is the nucleon mass. 

In Eq. ([3|, L{x) is a fiavor-indepcdent function, repre- 
senting the effect of the QCD interaction of the outgoing 
quark with the rest of the nucleon. The name "lens- 
ing function" has been proposed by Burkardt to denote 
L(x) [11) . Computations of the lensing function beyond 
the single-gluon approximation have been proposed in 
Ref. [12]. It is likely that in more complex models the 
above relation is not preserved, at least not as a simple 
product of x-dependent functions W. Nevertheless, it is 
useful and interesting to speculate on the consequences 
of this simple assumption. As a more refined picture of 
TMD and GPD emerges, it will be possible to improve 
the reliability of this assumption or eventually discard it. 
The present attempt should be considered as a "proof of 
concept" for further studies in this direction. 

The advantage of adopting the Ansatz of Eq. Q is 
twofold: first, it allows us to use the value of the anoma- 
lous magnetic moment to constrain the integral of the 
valence Sivers function; second, it allows us to obtain 
fiavor-decomposed information on the x-dependence of 
the GPD E and ultimately on the quark total angular 
momentum. This is an enticing example of how assum- 
ing model-inspired connections between GPD and TMD 
can lead to powerful outcomes. 

The Sivers function has been extracted from SIDIS 
measurements by three groups [T^16| . All of 

them assume a flavor-independent Gaussian transverse- 
momentum distribution of the involved TMD. Although 
this is an oversimplification, we adopt the same choice. 



At the starting scale Qo and following the notation of 
Ref. [19j . we use the unpolarized distribution and frag- 
mentation functions 



i^?(z,Pl;Qg)= ^°^^'^"^ -^l/<^^ 



n{Pl) 



(6) 



where z is the fraction of the energy of the fragment- 
ing parton a carried by the detected hadron. For fi{x) 
we use the MSTW08LO set [H], for D1{z) we use the 
DSS LO set ^8j. We fix the width of the transverse- 
momentum distributions for the initial parton and final 
hadron, respectively, as 



(ki) = 0.14 GeV^ (Pi) = 0.42 ^"•^'^(l - z 



\0.37 



GeV 



(7) 



These parameters have been implemented in the HER- 
MES gmc_trans Monte Carlo generator and are known 
to give a good description of the HERMES data [20] . In 
principle, these functions should be evolved according to 
TMD evolution j21j. However, we choose here to imple- 
ment only the evolution of their collinear part. 

Neglecting the contribution of heavier c, b, t flavors, we 
parametrize the Sivers function in the following way (in- 
spired by jT5]): 



/j^j. (x, K_|_; Qo) — /it l^'^Qo) 



Ml + (fcj) 
TiMl{kl) 



g-kl/Mt^-ki/{ki) 



(8) 



where Mi is a free parameter related to the width of the 
transverse-momentum distribution, and 



/it "{^'iQo) — ^ "v2e 
1 
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CV2e ,,f^!.. (1 - x) flix- Ql 
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Note that at Qo we establish a relation between the Sivers 
function for the combinations g^, q, and the correspond- 
ing unpolarized PDF, at variance with what has been 
done in the literature [15j[T6]. This will turn out to be im- 
portant when establishing a relation with the anomalous 
magnetic moment, since it guarantees that the valence 
Sivers function is integrable at any scale. We multiply 
the unpolarized PDF by (1 — a;) to respect the predicted 
high- a; behavior of the Sivers function [22]. We intro- 
duce the free parameter a"^" to allow for the presence of 
a node in the Sivers function at a: = a*^" , as suggested by 
diquark model calculations O [10] and phenomenological 
studies [13 (see the discussion in Ref. [23]). We imposed 
constraints on the parameters C" in order to respect the 



positivity bound for the Sivers function [2S], neglecting 
the contribution of the helicity distribution gi{x) (as in 
Ref. [H]). For the gluons, we assume the same func- 
tional dependence of the sea quarks, Eq. (10), with the 
replacement q ^ g. 

Also for /j^, we neglect the effect of TMD scale 

evolution [26]. We assume that f^j^ (a;;Q^) evolves 
in the same way as fi{x;Q^), based on the results of 
Refs. [27] [25] (note however that a slightly different re- 
sult has been obtained in Ref. [25]). 

In conclusion, we describe the SIDIS Sivers asymmetry 
in the following way: 
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and Pt is the modulus of the transverse momentum of 
the detected final hadron in the lab frame. 

For the lensing function we use the following Ansatz 



L{x) 



K 

{l-xY' 



(13) 



The choice of this form is guided by model calcula- 
tions [5HTU]. by the large-cc hmit of the GPD E [22], 
and by the phenomenological analysis of the GPD E pro- 
posed in Ref. [30] . We checked a posteriori that there is 
no violation of the positivity bound on the GPD £'*" as 
expressed in Ref. |31j . again neglecting the contribution 
of gi {x) . The nucleon anomalous magnetic moments are 
computed as 
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2£:"" (x, 0, 0) - E"^- {x, 0, 0) - E"- {x, 0, 0) 
2E'^- (x, 0, 0) - £;"" (a;, 0, 0) - E'" (a;, 0, 0) 



(14) 



We perform a combined x^ fit to 105 HERMES proton 
data [32], to 104 COMPASS dcuteron data [33], and to 
8 JLab neutron data |34j , of the Sivers asymmetry with 
identified hadrons. We sum the statistical and systematic 
errors in quadrature and neglect the experimental nor- 
malization uncertainty. Since the HERMES and COM- 
PASS data are presented as three projections of the same 
data set (binned in three different ways: in x, z, Ph±), 
we consider all three projections but we multiply their 
statistical errors by a factor v3 and we divide by 3 the 
number of these bins (105 and 104) when counting the 
number of degrees of freedom. The anomalous magnetic 
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1.591 ±0.009 


0.054 ±0.107 


-0.083 ±0.122 


Ml [GeV] 


K [GeV] 


?? 


a"" 


0.346 ±0.015 


1.888 ±0.009 


0.392 ± 0.040 


0.783 ±0.001 
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TABLE I: Best-fit values of the 8 free parameters for tfie case 
C" = C"" = 0. Tfie final xVdof is 1.323. Tfie errors are 
statistical and correspond to Ax^ = 1 



moments are known to a precision of 10~^ or higher |35j . 
However, given the typical uncertainties on PDF extrac- 
tions, our computation of k, is affected by a theoretical 
error of the order of 10"'^. Therefore, for our present pur- 
poses we take kP = 1.793 ± 0.001, k" = -1.913 ± 0.001. 

We started from considering 15 free parameters. They 
are C^, C"'",a'^", with q — u,d,s, the gluon coefficient 
C^, Ml, the lensing parameters K and rj, and the scales 
Qo and Q^. However, after some explorations, we made 
a common set of assumptions in all attempted fits. In 
all cases, we fixed a"*"'*" — (no nodes in the va- 
lence down and strange Sivers functions, as suggested in 
Refs. [ilinillSllMI)- We also set C^ = (the influence 
of the gluon Sivers function through evolution is anyway 
limited). Finally, all fits indicated that Qq = Ql ~ 1 
GeV was an acceptable choice. Therefore, the actual 
number of free parameters is at most 10. In this frame- 
work, we conclude that it is possible to give a simulta- 
neous description of the SIDIS data and of the nucleon 
anomalous magnetic moments assuming the relation in 
Eq. (|3|. 

We explored several scenarios characterized by dif- 
ferent choices of the parameters related to the strange 
quark. We considered fits with fixed C^ = 0, or with 
fixed C^" = 0, or with both parameters free (but con- 
strained within positivity limits), or with both fixed 
C^" = C* = 0. In all cases, we obtained very good values 
of x^ per degree of freedom (x^/dof) between 1.323 and 
1.347. All fits lead to a negative Sivers function for u„ 
and large and positive for d„, in agreement with previ- 
ous studies [T^HTB] and with some models [55H55] . The 
data are compatible with vanishing sea-quark contribu- 
tions (with large uncertainties). However, in the x range 
where data exist, large Sivers functions for u and d are 
excluded, as well as large and negative for s. The Sivers 
function for Sy is essentially unconstrained. The param- 
eter Ml is quite stable around 0.34 GeV, as well as the 
strength of the lensing function K around 1.86 GeV. The 
parameter rj is typically around 0.4 but can vary between 
0.03 and 2. The node a"" appears only above x ~ 0.78. 

We now discuss in detail the case with fixed C^" = 
C* — 0, because it gives the best x^/dof (1.323) and sug- 
gests that it is possible to fit the present SIDIS data for 
Sivers asymmetries in kaon emission without the strange 
contribution to the Sivers function. The best-fit values 
of the parameters are listed in Tab. |l] together with their 
statistical errors corresponding to Ax^ = 1. 




FIG. 1: The function xf-^j^ {x; Qo) (see text) as a function 
of X at the scale Qo ~ ^ GeV for a = u,d,u,d from top panel 
to bottom, respectively. The uncertainty bands are produced 
by the statistical errors on the fit parameters listed in Tab. [l] 



In Fig. IT] we show the corresponding outcome for 

^fiT K^'tQ'o) with a = u,d,u,d. The Sivers functions 
for s, s vanish identically. The uncertainty bands are pro- 
duced by propagation of the statistical errors of the fit 
parameters including their correlations, and correspond 
to Ax^ ~ 1. Our results are comparable with other ex- 
tractions of the Sivers function [T31 HH H^ . They are also 
qualitatively similar to the forward limit of the GPD E 
extracted from experiments ^U^ ISH EHl HO] • 

We can now compute the contribution to the anoma- 
lous magnetic moment of each valence quark flavor qy 
using Eqs. (14). We obtain 



= 1.673 ±0.003: 



-0.011 



-2.033 ±0m2toml, 
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The first symmetric error is statistical and comes again 
from the errors of the fit parameters (Ax^ = 1). The 
second asymmetric error is purely theoretical. It is com- 
puted by considering the other possible scenarios (cor- 
responding to different choices for C" and C"*) which 
give good x^ fits as well. However, a precise estimate of 
this error can be obtained only by performing a neural 
network fit [41] . The strange contribution to the anoma- 
lous magnetic moment is negligible, because the positiv- 
ity bounds severely limit the Sivers function for s and, in 



turn, also E'^^ and k*" . Our results are similar to other 
estimates of the strange Pauli form factor [42l [43] and 
lattice QCD calculations [iilliS] . 

Using Eq. (IT]), we can compute the total longitudinal 
angular momentum carried by each flavor q and q at our 
initial scale Qj^ = 1 GeV^. Using the standard evolution 
equations for the angular momentum (at leading order, 
with 3 flavors only, and Aqcd = 257 MeV), we obtain 
the following results at Q^ — 4 GeV'^: 



J" = 0.229 ± 0.002 
J'' = -0.007 ± 0.003 

J — U.UUO_Q QQg, 



+0.008 
0.012! 



J" = 0.015 ± 0.003 



+0.001 
^0.000' 



+0.020 
0.005' 



J^^^ 0.022 ±0.005to:S^oJ, 

/s - 006+°-°°° 
J — u.uuo_o 005. 



As before, the flrst symmetric error is statistical and re- 
lated to the errors on the flt parameters, while the sec- 
ond asymmetric error is theoretical and reflects the un- 
certainty introduced by the other possible scenarios. In 
the present approach, we cannot include the (probably 
large) systematic error due to the rigidity of the func- 
tional form in Eqs. (|8|)-(10), (13). The bias induced by 



the choice of the functional form may affect in particu- 
lar the determination of the sea quark angular momenta, 
since they are not directly constrained by the values of 
the nucleon anomalous magnetic moments. Our present 
estimates (at Q^ = 4 GeV^) agree well with other anal- 



yses 001 1311 [Ml EOl iSl iZl- It indicates a total contri- 
bution to the nucleon spin from quarks and antiquarks 
of 0.271 ± 0.007j:o:o28' of which 85% is carried by the up 
quark. 

In summary, we have presented a determination of 
the quark angular momentum assuming a connection be- 
tween the collinear limit of the generalized parton dis- 
tribution E and the Sivers transverse-momentum distri- 
bution. We have shown that it is possible to flt at the 
same time the nucleon anomalous magnetic moments and 
data for semi-inclusive single-spin asymmetries produced 
by the Sivers effect. Several different scenarios produce 
equally good x^ fits. Our strategy opens a plausible way 
to quantifying the quark angular momentum, and im- 
poses additional constraints on the Sivers function. 
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